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Abstract 

We present analytical results both in momentum and position space for the massless 
correlators of the vector and scalar currents to order as well as for the tensor 
currents to order af. The evolution equations for the correlators together with all 
relevant anomalous dimensions are discussed in detail. As an application we present 
explicit conversion formulas relating the MS-renormalized vector, scalar and tensor 
currents to their counterparts renormalized in the X-space renormalization scheme 
more appropriate for lattice calculations. 
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1 Introduction 



Correlators of gauge invariant quark currents are important objects in QCD. 
It is enough to mention that the correlator of two vector currents is directly 
related to the famous ratio R(s) = ^( e + e ~^^+ r / ° n )' > > t ne Adler function and the 
decay widths of the Z-boson and the r-lepton (for a review see, e.g. PQ). 

To be specific, let us consider a correlator: 

U(q) = i [ dxe^(Tj(x)j\0)), (1) 
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with j being a gauge invariant local operator. The polarization operator 11(g) 
satisfies the standard dispersion relationQ 



roo p(s) 

H(q) = / ds — subtractions, 

Jo s — q 2 



The subtractions on the right-hand side of Eq. (T5]) are necessary as they remove 
an additional divergence coming from the vicinity of the region x ~ in the re- 
integration in (pQ). The structure of the correlator (pQ) is significantly simplified 
if the momentum q is considered as large compared to (active) quark masses. 
Setting then all quark masses to zero one can describe the general structure 
of the correlator in pQCD as follows: 

n(g) = (Q 2 ) d - 2 f £ n^OO (i*^V- 

Here Q 2 = —q 2 , a s = 24 , \x stands for the renormalization scale and d is the 
(mass) dimension of the current j. 

In some applications it is useful to deal with the correlators in position space 
(see, e.g. p|3|r3|5|6] and below). Using a text-book formula for the massive 
scalar propagator in position space: 

1 



A(x,s) = / dq = zKi(z), z = V-x 2 s, (3) 

lyzTT) Jo s — q z —47T z x z 

(with K\ being a modified Bessel function) we arrive at a well-known repre- 
sentation for II: ^ 

n(a;) = / dsA(x,s)p(s). (4) 



o 

It should be stressed that the spectral density p(s) does not depend on the 
non-logarithmical contributions to the sum in ([T]) (that is those proportional 
to the coefficients 11^ with k = 0). Thus, the full correlator in position space 
considered as a function of x (defined for all x with x 2 ^ ) also does not 
depend on non-logarithmical contributions to 11(g). 

In general, the operator j is not scale-invariant (equivalently, has a non-zero 
anomalous dimension). The renormalization of the operator and the position 
space correlator look as follows^!] 

./ U(x) = Z 2 U (x), (5) 

where jo and n stand for the corresponding bare quantities. 



1 For simplicity we assume that the current j is a Lorentz scalar. 

2 Note that the momentum space correlator is renormalized in a more complicated 
way due to the UV divergence at small x. The corresponding formula is given below 
in Section [2j 
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An important feature of the position space correlators is that they can be 
directly computed non-perturbatively on the lattice by Monte Carlo simula- 
tions (see, e.g. [T|I5|9|TU] ) . Their long-distance behavior is governed by the non- 
perturbative features of the underlying field theory, QCD. On the other hand, 
due to asymptotic freedom, their short- distance behavior can be described 
by perturbation theory and operator product expansion (OPE). A meaning- 
ful comparison of perturbative results at short distances with their lattice 
counterparts requires, obviously, the use of one and the same renormalization 
prescription in the common case of scale-dependent operators. While minimal 
subtraction schemes (MS and its relatives [TO12II13] ) are certainly preferable 
for perturbative calculations, they, clearly, can not be implemented on lattice. 

A solution of the problem is based on the use of an intermediate renormaliza- 
tion scheme, with the renormalization conditions imposed directly on quark 
and gluon Green functions computed in a fixed gauge and for a particular 
configuration of external momenta |14j . 

A convenient intermediate scheme for the renormalization of the quark current 
operators has been developed in J9J. It is based on the study of the correspond- 
ing position space correlators and is called the X-space scheme. The conversion 
formulas between the MS and the X-space scheme have been elaborated in [9J 
to the next-to-leading order. 

Recently there has been a lot of progress in computing higher order correc- 
tions to the vector and scalar correlators within perturbative QCD both in 
the massless limit as well as for the general case of massive quarks. Both cor- 
relators are now known in momentum space to order [T5fl6|ll7|[T8] (the real 
and absorptive parts) and even, partially, to order [T9|20] (only the absorp- 
tive part in the massless limit). The situation is not so good for the tensor 
correlator, which is known completely to order a 2 s in the massless limit only 

In addition, the anomalous dimensions of the scalar and tensor currents are 
known to order a'l [22,23 24] (vector and axial- vector currents have identically 
vanishing anomalous dimension due to the corresponding Ward identity). 

The aims of the present paper are: 

• To compute the order contribution to the tensor correlator (only absorp- 
tive part in the massless limit). 

• To summarize available momentum space results for the scalar, vector and 
tensor (all massless) correlators and corresponding anomalous dimensions. 

• To discuss in detail the evolution equations for all three correlators. 

• To present full results for the correlators in position space, namely: scalar 
and vector to (and including) order a A s and tensor to (and including) order 
a? s computed within massless QCD. 
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• To construct the N 4 LO conversion formulas between MS and X-space renor- 
malized scalar and vector currents as well as the N 3 LO ones for the tensor 
current. 

• To study the stability of the conversion formulas with respect to higher 
order (not yet computed) perturbative corrections. 

The plan of the paper is as follows. In the next section we discuss our con- 
ventions and the definition of the X-space renormalization scheme as well as a 
version of the MS-scheme — the MS-one which seems to be more convenient 
for renormalization of the position space correlators. Sections [3] and H] list all 
available results for (massless) quark currents correlators in momentum and 
position space respectively. In Section we try to provide the reader with the 
concise bibliographical information about the origin of the results collected in 
the two previous Sections as well as about the main technical tools employed 
in the corresponding calculations. Conversion formulas between the X-scheme 
and MS/MS schemes are discussed in Section |6j In the last section [7] we sum- 
marize the content of the paper. 

In addition, there are three appendixes. In Appendix |A] we spell out the rules 
which we use to construct the Euclidean correlator from its Minkowskian 
counterpart. Appendix [B] provides the reader with necessary information on 
the Fourier transformation. Appendix O lists various anomalous dimensions 
relevant for the RG evolution of the quark current correlators. 



2 Quark current correlators in momentum and position space 

In this section we outline our conventions and recall the definition of the 
X-space scheme as presented in Ref. [9]. Our discussion will focus on the cor- 
relator of scalar currents first. The generalization to other Lorentz structures 
is straightforward, we will comment on it towards the end of the section. 

The scalar correlator in momentum space is defined as 

Ti S {q) = iJ dxe^(Tj(x)j(0)) (6) 

with j = ipttp. For space-like momenta we can express the correlator in terms 
of the Euclidean momentum Q. In what follows we will work exclusively with 
Euclidean correlators. Our procedure of obtaining Euclidean correlators from 
Minkowskian ones is described in Appendix [A] The correlator considered in 
position space reads 

U S (X) = (j(X)j(O)) (7) 

with a Euclidean separation X. We work in the chiral limit with = 0. 
Note that diagrams with purely gluonic cuts do not contribute to the scalar 
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correlator (in the assumed massless limit) 



2.1 Momentum space 



We denote the MS renormalized momentum space correlator at the scale \i by 

n s (Q, /i) = (z s ) 2 n s (Q) + z ss (^r e Q 2 , (8) 

where Uq(Q) is the bare scalar correlator. Note that in addition to the mul- 
tiplicative renormalization with Z s there is a subtractive counterterm Z . 
The corresponding renormalization group equation reads 



_d_ 

djl r 



with the anomalous dimensions 
2 d\ogZ s 



d/i : 



2 ' 
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A YSS 

ss (2 1 s + e)Z ss . 



d/i 2 



(9) 



(10) 



Using the solution to the renormalization group equation Q, we can evolve 
the correlator from one scale fio to a different scale \L\\ 



( a s {m) 



n 5 (Q,/ii) =exp 



\a a (Mo) 



dz 2 7 5 (z) 

~ W) 



A(/ii,/X ) 



a s (^o) 



dz j bb (z) 

~ W) 



t 



exp 



(n 5 (Q,/io) + Q 2 A(/i 1 ,/i 

dz' 2 1 s {z r 
~V f3(z') 



(11) 



\ a a (/xO) 



where a s = a s /ir = g 2 /(47r 2 ), g is the strong coupling constant and the (3- 
function /3(a s ) is defined as 



d/i 2 



J+2 



:i2i 



i>0 



While it is of course possible to recover logarithms explicitly using this solu- 
tion, it is more convenient to rewrite the renormalization group equation into 
a differential equation in I^q = \og([i 2 /Q 2 ) for this purpose: 



-^fI s (Q) = 2rTl*(Q) + ^Q 2 - a/n s (Q) 



(13) 



This equation can be used to iteratively reconstruct the logarithmic parts 
of II s (Q). Explicit formulas for the anomalous dimensions and the QCD (3 
function are given in Appendix O 
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2.2 Position space 



In principle, the discussion of the preceding paragraph can be directly trans- 
lated to the position space correlator. It is, however, convenient to use a mod- 
ification of the MS scheme that is a bit different from the traditional MS 
convention. The reason for this is that in the MS scheme logarithms in posi- 
tion space naturally appear in the forrrJjO 

\og(^-\+2 lE . (14) 



We can transform these to the simpler form \og{n 2 X 2 ) with a shift in the 
renormalization scale: 

H ->■ 2e" 7 > « 1.12/i. (15) 

The shifted [i defines a new modified MS scheme which we call MS. The 
relation between MS quantities and their MS counterparts is of course very 
simple: 

fl s (X, ii) = U s (x, 2e~ 7B /i) , 5 a (/i) = a s (W^/i) , (16) 

Using the evolution of the strong coupling constant we can also relate the MS 
coupling to the MS coupling at the same scale: 

a B (n) =a s (/i) j 1 - a B (n)lpo + a 2 s (fx)l((3 2 l - ft) + a 3 s (fx)l (-ft¥ + -ftftZ - ft) 



. i /:; ^3 „ 2 a j2 , o ( Pi 



PW - jPofal 2 + 3 I y + ftft 1 I ~ ft 



5/5^ 



+ 0{a b J 



(17) 

where I = 2 (log(2) — 7^). 

In position space there is no additional subtractive renormalization. Hence, 
the renormalization group evolution simplifies to 



z (5(z) 

K a s {no) 

The evolution equation for the MS scheme is obtained from Eq. ( fl8i) by simply 
replacing MS quantities by their MS counterparts. 



See Appendix IB. 31 for more details. 
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2.3 The X-space scheme 



The X-space renormalization scheme is defined by fixing the correlator of the 
normalized current jx = Z x j at a separation X to its value in the free 
continuum theory [9]: 

n|(x ) = (z x fu s (x ) = n s (x )| free . (19) 

This prescription can be readily implemented both in lattice and perturbative 
QCD. In perturbation theory the free theory value of a correlator is obviously 
just the leading order contribution. 



2.4 Other correlators 



In addition to scalar correlators, we also consider correlators of vector, ten- 
sor, pseudo-scalar and axial- vector quark currents. In position space these are 
defined as 

u p (x) =(j 5 (x) j5 (o)> , n*{x) =0,5(x)j, 5 (o)> (20) 

with 

js = #75^, jp = $1^, 3yh = ^7^75^, 3p.v = V><V^ ■ (21) 
Since we work with = 0, the results for the pseudo-scalar correlator will 
be the same as for the scalar correlator. 

Except for two small points, the entire discussion for the scalar case also holds 
for the more complicated Lorentz structures. First, in contrast to all other cor- 
relators, the vector and the axial-vector correlators do receive contributions 
from diagrams with purely gluonic cuts. We choose to neglect them in this 
work. This implies that also vector- and axial-vector correlators coincide. Sec- 
ond, it is not possible to naively renormalise the vector correlator according to 
the X-space condition (Eq. ([19]) ). The reason for this is that in position space 
its tensor structure varies between different orders of perturbation theory. We 
choose to renormalise the trace of the vector correlator instead. 



3 Momentum space correlators: results 



In the following two sections we present the results for the correlators both in 
momentum and position space. All results with their explicit renormalization 



7 



scale dependence can also be retrieved from 

http : //www-ttp .particle .uni-karlsruhe . de/Progdata/ttplO/ttplO-42/ 

We list the correlators in momentum space at the scale fi 2 = Q 2 , where all 
logarithms vanish. Results for arbitrary values of \x can be recovered by solving 
the renormalization group equation (i.e. by using Eq. ( ITT]) or ffT3]) ). Note that 
the anomalous dimensions listed in Appendix O allow the reconstruction of 
all logarithms at one order higher^ i.e. at order a\ for the vector and scalar 
correlators and at order o? s for the tensor correlator. 



n s W) = -^ 2 (i + £^ (n) ' s <) , 

C (1) ' S = ^-2Cs, 

17645 353 , 1 25 . / 511 2 \ 

215626549 1789009, 1639 , 2 1645, 73565 , 

248832 3456 s3 32 S3 1152 s4 1728 s5 

325 , 665, / 26364175 22769 , 5 , 2 53 

H Ce C7 + n/( 1 C3 - -Cl 

192 s 72 s ; V 373248 864 s 6 3 48 

1889 \ 2 /499069 157 1 5 \ 
+ "432" °7 + nf 1559872 " 1296 s3 + 48° " 18 C V ' 



5 / 00 

nJL(Q) = ^(-<? 2 V + QM ( 1 + £ c< " )- °' 

41927 829 . _ . / 3701 19 



( 3701 19 ,\ 

G + ^ + rifl 1 C3 

S3 S5 216Q 15 S3y 



C (2), 

1440 30 

,,, v 31431599 624799, 99 , 2 11, 349, 133, 

69120 1440 S3 2 S3 16 s4 32 S5 12 S7 
/ 1863319 174421 2 11 109 \ 

+ n/ 1 " ^4560~ + ^32T C3 " Cs " 48 C4 + T^ Cs ) 
2 / 196513 809 1 \ 

+ Hf V 155520 ~ 1080 s3 ~~ JT 5 / ' ^ ^ 
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& = (l + £ C7W'*a») , ^ = (l + £ D^aJ 

, 2U 556475 2657, 7 A 25, / 667 32 \ 

593 _ 12 
180 5 s 

, 2 w 566777 265 A 7 , . / 1333 64 \ 

£>(2),t _ a H /■ + 5 a + n 1 A ^4) 

19440 9 S 60 s S 'V 810 45 s J V ; 



with 



a 



T^Jpa ^pp^va ^pa^vp > 

T>(2) //ON _ QfJ-Qp f _ QpQcr r _ QuQp r QvQu r /nr\ 

1 pupaW) ~ Q2 °w Q2 ° V P Q2 ^ a ' Q2 W> 1 



and n/ active quark flavours. 



4 Position space correlators: results 



The position space results are obtained by four- dimensional Fourier transfor- 
mation (see Appendix [B]) of their momentum space counterparts. As it was dis- 
cussed in Section [1] the position space correlators are not sensitive (at x 2 ^ 0) 
to the constant (non-logarithmic) contributions to the momentum-space ones. 
Thus, the knowledge of the full vector and scalar momentum space correlators 
at order oP s plus the 0(ofy anomalous dimensions 7 55 and ^ vv from Appendix 
ICl allows us to present below the vector and scalar correlators in the position 
space at order a*. Similarly, the use of Eq. ([TBI and the 0(al) anomalous di- 
mensions 7W ,TT and ^ 2 ~)' TT computed by us (see eqs. in Appendix O) result 
to the full 0(al) results for the position space tensor correlator. 

We present the results in the MS scheme (see Eq. (fT6l) ) at the scale /i 2 = 
which correspond to MS results at the scale jj 2 = ^e" 27E . Results at an 
arbitrary scale fi can again be obtained with the use of the renormalization 
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group evolution (Eq. ([IS]) ). The correlators in position space read 



Tl 4 



£(2), 
£(3), 

£(4), 



n=l 



3' 

817 39 / 23 

144~T C3 + 7 H"72 
150353 5125 . 815 



5184 

2 



51 ^ + 12 ^ 5 + n/ V 3SSS 



13361 3 A 5 A 

+ J 6 " 5 C * 



+ n 



383 8 

1 C3 

11664 27 



+ 



22254833 592067 



458425 



497664 5184 
1375. 178045. 
^-Ce ^^C? + n f 



432 



,2 265, 

C 3 + T^Ci - 



607225 



864 



8775605 392129 



768 a " J V 373248 5184 
955 2 6731 45695 2875 665 

~ T6" C3 ~ ~576 U + ~21^ C5 + W Ce + TP 7 
,/ 224695 11263 , 5 ^ 25 , 6515 _ 



+ n 



2239488 1296 
3/ 6653 173 u I 



6 



96 



432 



f V 559872 1296 



144 C4 + 18 Cs 



7T J 

00 



C™ = 1 + C (n) ' v 



Y 



x 2 



~n 



£(2), 
£(3), 



n=l 
1, 

61 



n=0 



6 

7309 



iic 3 + M-| + fc 3 



989 275 / 2617 47 25 

48 6~ Cs + ~6" Cs + Uf { ~ 144" + Y Ca " Y Cs 



2 /277 8 A 

+ n, (3 

/ V648 

57640705 



27 

278401 , 5445 

Cs + — 5— CI 



133705 



20736 108 

/ 10278875 50705 , 

+ n f( ~ ~^nr + "77^3 



7315 



+ n 



\ 20736 144 
2 /1554751 10691 



288 ' 48 
2 65975 665 
55C 3 + ^-Cs + ^ t 



/ V 62208 
503 



5 2 _ 1315 

864 ^ 3 + £) 3 108 ^ 5 



+ « - 



2 _5_ 

1458 + 2T 3 + 18^ 5 
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D (0),v = D (l),v = Qj 

24 + 36 ' 

101 121 /587 11 



~ (3]v 101 121 /587 11. \ 2 / 1 l.\ 

~ (4) ,„ = _616333 mi _ 3025 /111409 _ 1607 275 

2304 4 s 48 ; V 2304 36 36 

2 / 54373 79 25 \ 3 / 325 5 \ 

+ M ~ 20736 + 36 Cs " 108 C V + "A 7776 " 162 C V ' 



C(i),t = 2 ; 

~ (2)t 5303 143 . / 443 2 A 

° (l =132-- l^ C3 + ^(- 648 + 3 C3 
~ (3)jt _ 8439437 77977 29 2395 

~ + 46656 48l) 3 + 54^ 4 + "l08^ 5 
/ 246515 1985 5 25 

+ Uf \ T166T + W C3 + 18 C4 " Y Cs 



with 



(27) 



2 /18287 4 \ 

+ "/(34992 - 9 C3 )' (28) 



, a,(i/V^) ~J_ 

7T f ( k 

k=l 

T (2) fy\_ X^X a X U X P X V X C 

and active quark flavours. 

In numerical form for rif = 3 the results read 

fl 5 (X) = — 7 (l + 0.67 S s - 16.3 5^ - 31 ^ + 497 a*) , (30) 

tt 4 (x 2 ) V ' 

KM = (% - ^r) (1 + a- - 2 - 5 a2 - 4 - 4 a2 + 66 

- V ( - 0-375 S 2 . + 0.63 a 3 s + 7.0 af) , (31) 

H^PO = 7^T3 fe^ - T£) CT (X)) (l + 2 5 S + 3.1 ~al + 6.6 a?) . 

(32) 
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This can be compared to the MS results at the same scale /i 2 = 1/X 2 : 

q 

U S (X) = 3 (l + 0.20 a s - 18.5 a 2 - 11 of + 579 a 4 ) , (33) 

tt 4 (x 2 ) V J 

K( X ) = f¥ - i 1 + a * ~ 3 -° a " ~ 2 - 4 a ' + 74 a 4 



7T 4 



(X^ V 2 A 2 

- 0.375 a 2 + 1.0 a 3 s + 6.4 a 4 ) , (34) 



= -"^7^3 - I&W) (l. + 2.15 a s + 3.3 a 2 + 6.4 a^ 



(35) 

where we use the abbreviation a s = a s (l/ \fX^j /tt. In both schemes the a 4 

coefficients of the vector and the scalar correlator are quite large. In the MS 
scheme we observe a slightly smaller a 4 contribution; furthermore the coeffi- 
cients are a bit more uniform between different orders. This indicates a better 
behavior of the perturbative series in this scheme. 

It is remarkable that the Lorentz structure of the tensor correlator remains the 
same at each order while it begins to vary at order a 2 s in the vector correlator 
result. This change can be easily understood by considering the transversality 
condition 

d,(j,(X)j u (0)) = 0. (36) 

Together with the appearance of terms which are logarithmic in X also the 
Lorentz structure has to change so that the vector correlator remains transver- 
sal (see Eq. (1B.4|) in Appendix [B]) . 



As a consequence of this, it is not possible to renormalise the vector correlator 
in the X-space scheme according to the prescription Eq. (fl~9l . There are many 
possible generalizations leading to a well-defined renormalization prescription 
for the vector correlator. Among them, we choose to renormalise the trace 
of the correlator. For the other correlators the renormalization is straightfor- 
ward. We refrain from presenting the somewhat lengthy explicit results here. 
They can be easily constructed from the MS results (Eq. f[26]) - (T28|) ) and the 
conversion formulas listed in Section [6] (Eq. fj40|) - fj42|) ). 



5 Bibliographical and technical comments 



We start with some generic notes. First, the non- logarithmic part of a quark 
current correlator in momentum space is not physical within QCD as it re- 
quires additional UV subtractions beyond the ones associated to the coupling 
constant renormalization. 
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Still, the constant part is important because it provides us with a very conve- 
nient way to compute the Q-dependent (that is physical) contribution. Indeed, 
suppose that we want to get the Q-dependent contribution at order a™ to, say, 
the scalar correlator (the discussion below is general and applicable to every 
massless correlator). A direct calculation would require to deal with (n+1)- 
loop massless propagator- like diagrams contributing to II 5 . A better waj0 is 
to use the evolution equation ffT3"l) . Indeed, after a (trivial) integration of the 
right-hand side of the equation with respect to I^q we arrive at the conclusion 
that the Q-dependent part of II s is completely determined by the knowledge 
of three ingredients: 

1. the very correlator U s (including its constant, Q -independent part) and the 
beta-function to order a™ (both are contributed by n-loop diagrams); 

2. the anomalous dimension / -f s to order a™ ( n-loop diagrams); 

3. the anomalous dimension r y ss to order ((n + l)-loop diagrams). 

Moreover, it is a well-known fact that the methods of Infrared Rearrangements 
(IRR) and .Reoperation [3U] allow to reduce the problem of evaluation 
of a (n+1) loop contribution to an (arbitrary) anomalous dimension to the 
computation of some properly constructed set of n-loop massless propagator- 
like diagrams. 

Up to and including three loops the massless propagator-like diagrams can be 
computed easily with the FORM J3TJ package MINCER [52153] . The package 
implements the algorithm developed in [33] and is based on the use of the 
traditional method of integration by parts. Another powerful approach - 
the method of the Gegenbauer Polynomials in x-space (GPTx) [13J - - is less 
automated, but, sometimes, is applicable in cases with loop number exceeding 
three (see, e.g. [35J for recent spectacular examples). 

The only systematical way to compute massless propagators at four loops is 
based on the so-called 1 / D-expansion elaborated in [5oTI5T] and on the use of 
a special parametric representation of Feynman integrals [38.39 40J. Here by 
computation we mean the reduction to the corresponding master integrals, 
the latter have been computed analytically [JT] and numerically [42] . 

In Tabled] we display the bibliographic information about results listed in Sec- 
tions ISHl and in Appendix B. In addition we also mention the main theoretical 
tools used in obtaining these results. 

We want also to stress, that literally all 0(af) calculations listed in the Table 



4 The observation below was explicitly made (for the particular case of the vector 
correlator) in |25j . 
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a 2 s 




<4 


TT 5 

11 


Il , 7l91l MTNfFR 
ii/ p^jii. iviiiN v^njrv 


1 (_J /r~\ ill} nvn o n ci i~\n 
1 IjpZrUI , 1 / LJ C-A-JJclllolUll 






[27|T7], MINCER 


[ITJ, MINCER, IRR 


PPB], 1/D expan- 

cinn 

OlUll, 

-Reoperation 


Tm T\ S 


[971 MTNC!F,R 

IZj 1 1 . 1V-L-L1 M 


fTTl MTNC!F,R 

11 1 1 , 1 VI 111 V 1 ^11. 

i?*-operation 


1 Q \ I 1 I pvn anci on 

1 -L Cfl ; ± / J-^ C-iVJJdllOlvJll, 

-Reoperation 


11 


llt_)|lUpZfl|, 1V111N V^JjjJA, 


1901981 1 / D pvnansinn 






[25], MINCER 


[15|16j, MINCER, IRR 


|20|28], 1/D expan- 
sion, 

i?*-operation 


Tm T\ v 


[9^1 TRR 


rm mtnpfr 

i?*-operation 


) >4 T / / J ovr\ o n c i /~\n 
IZOI , 1 / LJ cXpdilolUll, 

i?*-operation 


11 


urn mtnpptj 






TT 


[21], MINCER 


present work, MINCER, 
-Reoperation 




imn T 


[ZI], MINCER 


present work, MINCER, 
-Reoperation 





Table 1 

Calculations of massless quark correlators in QCD: references and main theoretical 
tools. 



would be not possible to perform without a heavy use of the parallel versions 
of FORM, PARFORM [35133135] and TFORM [36]. Last, but not least, the 
generation of thousands of four-and five-loop input QCD diagrams have been 
conveniently done with the help of the (FORTRAN) program QGRAF [37]. 



6 Conversion between X-space scheme and MS or MS 

The natural scale for a transition between the MS and the X-space scheme is 

= 4^ ■ (37) 
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This is obviously equivalent to a transition between MS and X-space scheme 
at a scale /1q = -^-e -275 . Imposing the X-space renormalization condition 
(Eq. (USD) on 



7& \ 2 



n d (x ,/io) 



7& \2 



U 5 (Mq) 



n 5 (x ,/i ) 



(3* 



with 5 e {V, S, T} directly yields the desired ratios jH] 



Z 5 (fi ) Z s (jl 



Z' x 



7 s 



U s (X ,fl ) 



\ n 5 (x 



free 



(39) 



between the renormalization constants in the different schemes. We obtain 



z s {h) 



Z x 



l + £^)> s a s (/V) r 



71=1 



6 (i), 



3' 

89 39, / 23 

G + n f [ 

32 4 s ; V 144 

140741 2387 815 

10368 ~ "l^P 3 + ~24^ 5 + Uf 
2 / 383 4 , 

+ ni 1 C3 

J V 23328 27 s 



12947 19, 

+ ™C 3 



7776 72 



13901515 4393 



995328 
1 3 75 

<6 



288 
178045 



208679 



12 



^ 25 > 
C4 - tt:Ce 



432 



64 
851 

32 S3 



(r + n 



265 

<3 + "So^ 4 ~ 
8029687 



626785 

1728 
418799 



C5 



1536 a ' J V 746496 10368 
6571 45895 2875 665 

~ 1152 C4 + ~432~ C5 + ^76~ C6 + 144 Cr 
257315 11273 13 2 25 6515 

4478976 + 2592 ^ 3 + 36 s3 + 192^ 4 ~ 864 ^ 5 
0/ 6653 173 , 1 , 5 > 

+ n 3 f C3 H C4 H C5 

7 V 1119744 2592 s 288 36 



(40) 
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^M = (^r i = i + £^s s ( / 2 r ) 

L X n=l 
Jf(l),« - I 

2' 

r / 2 s„ 97 XX / X 1 \ 

,„ = 14881 _ 119 275 / _ 5395 47 _ 25 , 

288 2 s3 12 S5 1 V 864 9 s X8 S V 
2 / 181 2 \ 

+ n Ax296 "27 C3 J' 
x(4),« _ 34042561 29437X 5203 2 _ 2X2905 _ 73X5 

d " 4X472 432 ^ 3 + ~X() 3 576 ^ 5 ~ ~96~" 

/ 6096767 78X9, 77 , 2 79775, 665 , \ 

+ n f [ 1 (3 H C5 H C7 

; V 4X472 96 s 3 S3 864 s5 X44 s 7 

2 /889699 2899, X3 >9 1415 , \ 

+ m C3 H Cl C5 

/ V1244X6 X728 36 3 216 / 

3 / X037 2 , 5,\ 

+ "H-Tl664 - 81 6+ 36 C5 )' (41) 



!^ = i + f; *(»>.< s , (Ao) », 

<y«'* = x , 

r / 2 w 487X X43, / 443 X, 

§m = * + _| — a 

864 36 ; V X296 3 / 

rf ow 79X3369 X8529, 29 , 2395, 

933X2 243 X08 2X6 



+ n f 



23854X 19X3 ^ 5, 25 , \ 9 /X8287 2 \ 

1 C3 H C4 C5 +ni C3 • 

23328 2X6 36 18 V / V69984 9 / 



(42) 



The conversion can of course be performed for an arbitrary scale fi with the 
help of the renormalization group evolution (Eq. ( 1X8]) ) . 

The new higher order results reduce the theoretical error significantly com- 
pared to the old NLO conversion formulas of Ref. [9]. As an example we show 
the values of Z 5 {2 GeV)/Z 5 x (X%) with 5 G {V, S, T}, X 2 = (1.5 GeV)" 2 at 
different orders in perturbation theory. We estimate the theory error from 
higher orders by performing the transition between the two schemes at some 
varying intermediate scale /x and evolving the result to the final scale of 2 GeV: 



Z s {2GeV) 



exp 



/a s (2GeV) 

' ' dzj s {z) Z\ii) 



z (3(z) I Z X (X^ 



(43) 



For the evolution, we use anomalous dimensions at one order higher than 
Z 5 {2 QbV) / Z x {Xq) , up to the highest available order of a 4 . In Figured] the 



X6 



scalar correlator 




N 
> 
N 




2 2.2 2.4 2.6 2.i 
H [GeV] 



vector correlator 




N 
N 




2 2.2 2.4 2.6 2.8 3 
H [GeV] 



tensor correlator 





2 2.2 2.4 2.6 2.i 
H [GeV] 



Fig. 1. Values of Z s (2 GeY) / Z & x (Xq) for vector, scalar and tensor currents with 
Xq = (1.5 GeV) -2 and nj = 3 quark flavours at different orders in perturbation 
theory. The error bars on the left-hand side are obtained by running the ratio of 
renormalization constants from an intermediate scale ^ with 1/Vq < [i 2 < 4/Xq to 
2 GeV. On the right-hand side, the dependence on the choice of [i is shown. 
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scalar correlator 




N 
N 



order 



order 




1.6 1.8 2 2.2 2.4 2.6 2.8 3 
\i [GeV] 



vector correlator 




N 
N 




tensor correlator 




N 
N 



2 2.2 2.4 2.6 2.! 
H [GeV] 




2 2.2 2.4 2.6 2.8 3 
H [GeV] 



Fig. 2. Values of Z s (2 GeY) / Z 6 x {Xq) for vector, scalar and tensor currents with 
Xq = (1.5 GeV) -2 and nj = quark flavours at different orders in perturbation 
theory. 
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values for this transition factor are plotted for a varying intermediate scale 
fi 2 between l/Xfi and 4/Xq and rif = 3 quark flavours. Figure [2] shows the 
same ratios of renormalization constants for rif = 0. As expected from the 
numerical formulas Eqs. - (JHSJ) the vector and scalar correlators receive 
large contributions at order a 2 s and a* for rif = 3. 



7 Conclusion 



In this work we have presented presently available information on three basic 
quark currents correlators — the scalar, vector and the tensor ones — con- 
sidered within the massless QCD. The correlators and the corresponding RG 
evolution equations have been studied both in the momentum and position 
space. Explicit conversion formulas relating the MS renormalized vector, scalar 
and tensor currents to their counterparts renormalized in the X-space renor- 
malization scheme are constructed. It is demonstrated that the new higher 
order results reduce the theoretical error significantly compared to the old 
NLO conversion formulas of Ref. [9]. 
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A Euclidean correlators from Minkowskian ones 

Let us consider a generic correlator defined originally in the Minkowskian 
space: 
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IW,..^,..(g) =i f M dx e^(Tj ulU2 Jx) jI, 2 .M) (A.l) 

= E^ 1 ,,..(^ 1 ,,..(g)n fe (-g 2 ), (A.2) 
k 

where the tensors ViU2 are made from the metric tensor g^, the vector 
q and the indexes //1//2 • • 

The corresponding Euclidean correlator in the momentum space is constructed 
as follows: 

rw 2 ...^ 2 ...(Q) = E T ^ 2 ...^,..(Q) MQ 2 ), ( A -3) 

with T^ ifl2 UiU2 (Q) being made from T^ W2 vxV2 (q) with the help of the re- 
placements q^ -Q^q Vi ->■ -Q Vi ,g^ -»■ -<W 

At last, the Euclidean correlator in the position space is de/in erf with the help 
of the Fourier transformation, viz. 

IW..™...(Q) = / dX e** (Tj nn jX) jl W2 .M) (A.4) 



B Fourier transformation in d dimensions 



Perturbative calculations are in most cases quite cumbersome in position 
space. A more convenient alternative is to perform the calculation in mo- 
mentum space, where exploiting invariance under translations (implying mo- 
mentum conservation) leads to great simplifications. In the end, the result in 
position space can be recovered by means of Fourier transformation. 



B. 1 Fourier transformation of bare functions 



In the current work, all bare Green functions have a power-like dependence 
on the momentum. For the transformation from euclidean momentum space 
to euclidean position space we use the following formula: 



FT ( (0 2 )~ r } — f ^ ^ - 1 T{i i- r) ( X2 \ ' (BY) 

It is straightforward to derive this formula by using Schwinger parametrization 
and Gauss integration. 
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For Green functions with a non-trivial Lorentz structure it is convenient to 
use momentum operators, i.e. 



->■ - 1 



d 



(B.2) 



FT (Q,Q u (Q 2 )- (r+1) 



(47T) 



u flV 

~2r 



2r J X 2 



(B.3) 



Note that forming the trace of the right-hand side of Eq. flB.3j) again yields 
the simpler result from Eq. (IB. II) . 



We also see explicitly that the Fourier transformation preserves transversality. 
Combining Eqs. ( IB.lj) and (1B.3j) we obtain 



FT 



QuQ 



d/j, FT 



Q 2 



QuQv 



(Q 2 



1 r(f -r) fx 



(4vr; 



I r(r) 



(B.4) 



x 



1 - — )6. 



2r 



J flU 



d/\ X^Xy 

~ 2r X 2 



Q 2 



(Q 2 



--d u FT 



Q 2 



(Q 



0. (B.5) 



B.2 Fourier transformation in four dimensions 



The Fourier transformation for the correlators we consider is regular at d — 4. 
This means that there are two ways of obtaining the renormalized correlators 
in position space. On the one hand, we can transform the bare Green functions 
in d dimensions and perform the renormalization afterwards. On the other 
hand, we can start from the renormalized correlator in momentum space and 
transform it in four dimensions. In the latter case we have to transform terms 
which are logarithmic in Q 2 in addition to simple powers of Q 2 . Rewriting the 
logarithms with the help of derivatives leads to 



FT 



(Q 2 ) r log s (Q'< 




FT 



r+S 



(B.6) 



In Table IB. II the resulting expressions are listed for all powers of Q 2 and 
log (Q 2 ) which appear in the momentum space correlators up to order a^. 
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Momentum space Position space 



log^/Q 2 ) 
logV/Q 2 ) 
logV/Q 2 ) 




logV/Q 2 ) 






iog 5 (/i 2 /g 2 ) 


_^ ( _16C 3 + 16Z M xC3-^ x 




Q 2 log(^ 2 /Q 2 ) 


8 

7T 2 (X 2 )3 




Q 2 logV/Q 2 ) 






Q 2 logV/Q 2 ) 






Q 2 logV/Q 2 ) 


_^ f ( 6 _ 8C 3_24/ mX + 15/ 2 x 


- Kx) 


Q 2 logV/Q 2 ) 


i ^ F (40C 3 + 12Z / ,x-16Z M xC3 


- + ItMjx - #x 



Fourier transformation of logarithmic terms for various powers of Q 2 and logQ 2 . 
Terms that do not contain logarithms vanish for non- negative integer powers of Q 2 . 
We use the abbreviation l^x = log(^ 2 ^ 2 /4) + 2^e- 

B.3 Logarithmic structure 

One-scale problems, like the correlators considered in this work, exhibit a very 
special structure of their logarithmic terms both in momentum and in position 
space. 

In momentum space in the MS scheme, logarithms have the form \og(fi 2 / Q 2 ) . 
They originate from terms of the form 



where I denotes the number of loops in the corresponding diagram. In order to 
obtain the logarithmic structure in position space we can again extract terms 
of the form y le from the Fourier transform of the left-hand side of Eq. ( IB. 71) . 
Expressing the Gamma functions in the Fourier transform in terms of expo- 
nential functions and polynomials in e we find 



Table B.l 




(B.7) 



( 



2 X 2 



) 



lc 




2 X 2 



) 



4 



1+le log 



4 



+ 2 lE +0(e 2 ), (B.8) 



for the general structure of position space logarithms. 
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C Anomalous dimensions 



In our conventions the momentum space renormalization group equations of 
the scalar, vector and tensor correlators read 



/i 2 ^n*(Q, /i) = 2 1 S U S (Q, p) + 7 ^Q 2 , (C.l) 
d 



/i 2 — n^(Q,/i) = 2 1 v Ul u (Q,f,)+ 1 vv (-Q 2 5 IIU + Q fl Q l/ ), (C.2) 

^^(Q'ri = «,(^) + (7 (1) ' tt t« ct + i^ tt t$ p aq))q 2 

(C.3) 



with the tensor structures 



T^Jpcr dfipfiva <^p.o^vp i 

T (2) (n \ _ QyQp x _ QtiQc r _ QvQp c , QuQa r (c 



The corresponding evolution equations in position space are analogous, but 
contain no subtractive anomalous dimensions. The anomalous dimensions are 
given by 



7>s)=0, (C.5) 

±U/I n=0 

7^ = 4, 

7r=4, 



125 11 



72 = 72--18 n " 
Vl/ 10487 110, / 707 110 A \ 77 , 

73 =^- + -9- C3 + n/ l _ 2i6 "17 Cs J - 972^ ' 
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74 ~ 41472 + 864 ^ 3 ~~ To ~ 288 ^ 5 

( 11785 58625 715 13325 

+ n/ l"~648 86T Cs + l8" C4 + ^32" Cs 

2 / 4729 3163 . 55 , \ , / 107 

+ n 2 f [ + Cs C4 + nil + Cs , (C.6) 

f \ 31104 1296 s 72 s / 7 V 15552 108 s ' ' 
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Additionally, the four loop QCD /3 function is required for the renormalisation 
group evolution (Eqs. (ITT]) . ( ITS]) ). In our convention it reads 



(5(a s ) = - Pna 
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The (3 function and the mass anomalous dimension 7 m were computed at four 
loop order in Refs. pE9l22ll23] . 
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